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Abstract 

An analysis has been made of the present situation with the high energy hadron- 
nucleus and nucleus-nucleus interaction models. As is already known there are 
inconsistencies in the interpretation of experimental data on the primary mass com¬ 
position, which appear when different EAS components are used for the analyses, 
even for the same experiment. In the absence of obvious experimental defects, there 
is a clear! need for an improvement to the exising models; we argue that the most 
promising way is to introduce two effects which should be present in nucleus-nucleus 
collisions and have not been allowed for before. These are: a few percent energy 
transfer into the EAS electromagnetic component due to electron-positron pair pro¬ 
duction or electromagnetic radiation of quark-gluon plasma and a small slow-down 
of the cascading process in its initial stages associated with the extended lifetime of 
excited nuclear fragments. The latter process displaces the shower maximum deeper 
into the atmosphere. 

1. Introduction 

Among the most popular interaction models which are used for the analysis of exper¬ 
imental data in the PeV region and beyond are QGSJET, SYBILL, NEXUS, DPMJET, 
VENUS and HDPM, coupled with CORSIKA and partly with the AIRES simulation 
codes ( see Bi and references therein for a description of the models ). Early versions 
of these models differed signihcantly in their description of the development of the atmo¬ 
spheric cascades. One of the possible ways to evaluate the quality of a model is to derive 
the primary mass composition from the analysis of different observables: a good model 
should obviously give the same composition for any set of observables. We used this ap¬ 
proach to interpret the data on the maximum development depth of the electromagnetic 
cascades X^ax and on the ratio of muon to electron sizes N^/Ne in terms of the mean 
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logarithm of the primary mass (InA)] the best consistency was obtained for the QGSJET 
model [^. A later analysis of the EAS hadron component conhrmed this conclnsion Q. 

Since then the models have been improved and the differences between their predictions 
rednced [0 . The precision of experimental data has been also improved and this offers the 
hope of fnrther testing the models and of determining the actnal valne of (InA). Some of 
the new measnrements in which Cherenkov light in the atmosphere is detected have given 
a deeper X^ax and hence a lighter mass composition than before [^, The extensive 
set of measnrements on the gronnd and relevant simulations made by the KASCADE 
collaboration indicated that there is a difference in the estimate of (InA) derived from 
different sets of observables 0]. At lower energies it was also difficult to get a consistent 
explanation of the measured muon and hadron trigger rates [§. All these facts provide 
evidence that further improvements to the models are urgently needed. 

Meanwhile, a few calls for a quite radical change of the interaction model have ap¬ 
peared, inspired mostly by the existence of the well known knee in the primary cosmic ray 
( CR ) spectrum at PeV energies. We discuss hrst the need for such a radical change and 
follow this with our own suggestion for modihcations to the model for a better description 
of the modern experimental data ( i.e. KASCADE, DICE, BLANCA et al. ). 

2. The possibility of a radical change in the interac¬ 
tion model 

2.1. The possibilities and general remarks 

A number of possibilities arise, and they will be considered in turn. Essentially they are 
of two types. The hrst considers that the primary spectrum ( represented in the usual way 
by a power law ) has an energy-independent exponent valid both below and above the PeV 
region where there is the knee in the spectrum inferred from the detected particles ’low’ in 
the atmosphere ( electrons, muons, hadrons ). The change in the interaction mechanism 
is then sudden ( at ~ 3 PeV ) for CR-air nucleus interactions. The second is of the same 
philosophy i.e. that the original production spectrum has a unique exponent, but the 
interactions are in the interstellar medium ( ISM ) or even in the intergalactic medium 
( ICM ) and not with air nuclei so that the spectrum of CR incident on the atmosphere 
already has the knee in it; there is no further radical change in the interaction mechanism 
for CR - air nucleus collisions. 
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The general argument against all these ideas is based on the observed anisotropy of the 
arrival directions of cosmic rays. The measured amplitude and phase of the anisotropy 
change in the vicinity of the knee and this gives support for its astrophysical origin 
i 0 Q, i-e. that the spectral shape arises because several sources, or types of source, 


contribute. 

Of relevance, too, is the cosmological aspect, viz the effect of a proposed radical 
interaction change in the early Universe when the average energy of all ’particles’ was 
in the PeV region. As we have pointed out already there are sound cosmological 
arguments against a radical change for P-P reactions in the PeV region. Some 10“^^ sec 
after the Big Bang the average temperature of the primordial plasma was about 10^®K 
and the mean energy per ’particle’ was of the order of a TeV in the centre of mass system, 
which corresponds to PeV energies in the laboratory system. If a radical change of the 
interaction occured at this energy changes to the subsequent evolution of the Universe 
would have taken place |^. Now the contemporary model (i.e. without a radical change ) 
fits the observed data on the relative abundance of the light elements and it is difficult to 
imagine this being the case were there to have been a radical change of the interaction at 
10-1^ s. 

Despite ’the case against’ the possibility of such a radical change it would have such 
important implications for particle physics that we make a detailed examination of some 
of the recent radical interaction models put forward to explain the knee. 

2.2. A radical change for the CR - air nucleus interaction 


Nikolsky has argued extensively that there is such a radical change. The latest idea 


energy of about 3-5 PeV the multiplicity of secondary particles increases dramatically. 
The attenuation of the cascade speeds up and the intensity of showers in the lower half 
of the atmosphere falls, thereby forming the knee. There is actually no detailed model 
to be checked, because no numerical estimates have been given of the basic interaction 
parameters and their energy dependence, therefore it is difficult to assess the idea quan¬ 
titatively. 

One remark should be made, however. There is a difference between the previous and 
present pictures of the interaction advocated by the author. The 1995 model had the need 


1311 is close to that of his paper of 1995 [H. The author assumes that at a primary 
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to withdraw a part of the cascade energy and to put it into some unobservable components 
( high energy muons, neutrinos etc. ) in order to form the knee. Now, in 2001, on the 
other hand, more energy is transferred to low energy particles of the hadronic cascade in 
the initial stages of its development. In our mind this is a backward step. The suggested 
increase of the energy transferred into the hadronic cascade in the atmosphere should be 
seen in the spectra of muons and Cherenkov light emitted by the charged particles in the 
atmosphere as an ’inverse knee or ankle’, i.e. as a decrease of the spectral slope. Nothing 
like that is actually observed. EAS muon size spectra at mountains: Tien-Shan |]^ p!6| , 

I and MSU 1^ all show 


EAS-TOP 


and at sea level: Ohya |T3, KASCADE 


the ’ordinary’ knee, i.e. an increase of the spectral slope. To the existing 5 Cherenkov 
light spectra: TUNKA H, HEGRA H, DICE §, BLANCA |[ and CACTI [H has 
been added another one which is based on the combined use of Cherenkov light and EAS 
electron size measurement [^]. All 6 demonstrate clearly the existence of the knee just 
like that in the particle size spectrum with a sharp increase of slope, which rules out 
the possibility of radical increasing the energy fraction deposited by the cascade in the 
atmosphere. 

The possibility of having no knee in the primary spectrum and hiding the missing 
energy in high energy muons and neutrinos ( an apparently attractive idea ) has also been 
discussed by Petrukhin [^. The proposed model suggests the production at threshold 
of one or a few hypothetical heavy particles in each interaction, which then decay into 
muons and neutrinos. The energy flux contained in his muon spectra is less than the 
missing energy flux, needed to form the knee, by a factor of about 20 even for the case of 
the maximum flux of high energy muons. This means that the bulk of the missing energy 
should be carried away by neutrinos. The author succeeds in making the energy spectrum 
of muons from the hypothetical decay much harder than that of the background muons 
from 7t,K decays and prompt muons. In this way it is possible to make the missing energy 
visible only as a very low flux of extremely high energy muons. However, such a model 
of the ’new physics’ ignores everything we know about the multiple production at high 
energies ( small cross-section for the production of heavy particles, conservation of lepton 
number in weak decays, etc.) and seems hardly credible. 


Another model of ’new physics’ discussed by Kazanas and Nicolaidis [28|, in spite of 
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some theoretical arguments, is actually purely phenomenological. In it one opens a new 
channel for multiple production and puts part of the primary energy into some unobserved 
particles (technihadrons, gravitons ), then, if one hts the energy dependence of the cross- 
section for the new channel and the energy fraction contained in the new particles, one 
can get a reasonable £t to the observed primary cosmic ray energy spectrum, including 
the knee. It is difficult to argue against the £t, but there seems to us to be no support 
for such speculation. 

2.3. Radical changes in the interactions in the interstellar 
medium 


It is, in principle, possible that the knee is due to interactions of some form beyond the 
atmosphere - in the ISM or IGM. If, as might be assumed by analogy with the ’X-particles’ 
produced by cosmic strings, the electromagnetic component is a prominent end-product 
of this interaction, then this hypothesis can be ruled out because the ensuing cascade in 
the Universe ( inverse compton scattering on the cosmic microwave background ) would 
produce a gamma ray intensity in excess of that measured ( see, e.g. for the related 
cosmic string arguments ). However, there are other possibilities for the radical change 
of the interaction, which can now be considered. 

The possibility of hnite neutrino mass has given birth to models in which the knee is 
caused by the interaction of the primary protons with massive relic neutrinos. Assuming 
that the mass of the electron antineutrino is 0.4 eV/c^, Wigmans converts some of the 
primary protons by inverse /3-decay reaction into neutrons and positrons, starting from 
an energy of about 4 PeV. Indeed, protons will disappear from the primary flux, but we 
ourselves argue that even if the shower from a secondary positron misses the detector, 
the shower produced by the neutron will be practically indistinguishable from a proton- 
initiated shower, because neutrons preserve more than 99% of the primary proton energy; 
thus EAS measurements would not see a knee. Taking into account the problem of the 
very small cross-section for the neutrino induced reactions, which requires us to assume a 
very dense halo of neutrinos gravitationally attracted to the cosmic ray source region ( be 
it the Galaxy or source within the Galaxy ), one should remark that the PeV neutron will 
decay again into a proton and a positron within tens of parsecs and the model loses its 
value. 
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Another model of the intermediate type was proposed by Ehrlich He used some 
theoretical ideas that the electron neutrino is a tachyon, which in fact does not contradict 
the existence of a hnite neutrino mass. In this case, starting from the same energy of 
4 PeV, the primary proton has the possibility to decay into a neutron, a positron and 
an electron neutrino. This decay might create a narrow peak of neutrons in the primary 
cosmic ray flux in the region close to 4.5 PeV, which can perhaps imitate the knee. All the 
scepticism put forward in the previous paragraph is applicable also in this case, however. 

The common dehciency of all the models on the radical change of the interaction is 
that no attempts have been made to reconcile it with the multitude of other experimental 
data, different from those used for the motivation of the model. The conclusion which 
can be drawn from this survey is that at the present time there is neither a reasonable 
model nor solid argument which justifies a radical change of the interaction. We conclude 
that such a change as is necessary for the interaction model should not be radical and 
not responsible for the origin of the knee, but only responsible for restoring consistency 
for all cosmic ray measurables. 

3. What kind of the modification of the interaction 
model do we need ? 


The variety of observed EAS components and precise measurements of their character¬ 
istics by the KASCADE experiment allowed this collaboration to derive the primary mass 
composition using the multivariate analysis of data. Remarkably, as has been pointed out 
already, it was found that not only do different methods applied to the same set of observ¬ 
ables give different results, but that there is a systematic difference between the results 
obtained using the same method, but different components. In Figure 1 the mean loga¬ 
rithm of primary mass (InA) derived from different sets of observables is shown vs. the 
so called truncated muon number Njf ( the hgure is taken from [0 ). Njf is the number 
of low energy muons ( > 0.25GeV ) contained between 40 and 200 m from the shower 
axis and is an observable adopted by the KASCADE collaboration as a measure of the 
primary energy, independent of the primary mass. They also use the number of higher 
energy muons ( > 2.4GeV ) Nf as well as some observables related to the hadron compo¬ 
nent for the analysis of the mass composition The fundamental problem is, of course, 
that all the (InA) values should be the same, at the same value of Njf, and they are not. 
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Figure 1: ( a,b, taken from 0 ). Mean logarithmic mass {InA) from the analysis of different sets of 

observables vs. IgNj^ ( QGSJET prediction ). The two dashed lines indicate alternative predictions of 
the Single Source Model with the upper line being our latest and ’best estimate’. The sets displayed 
in (b) do not include the observable Ng. It is seen that omitting Ng results in a heavier mass composition, 
(c) The ratio of mean values of {InA) from sets (b) and (a). We argue, in the text, that the ’true’ value 
of (InA) is a little below the crosses in Figure lb. 
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A number of features are of interest and, presumably, of some importance. 

(i) The dispersion of the individual {InA) values about smooth lines through them 
indicates that non-systematic errors are ~ ±0.05. 

(ii) The use of the electron size results in an average lighter composition ( Fig.la ). 
On the other hand, omitting electrons and using just muons and hadrons results in a 
heavier composition ( Fig. lb ). 

(iii) Considering the mean values from the two sets of data ( £lled-in circles ), the 
ratio of {{InA)) for (b) to (a) increases smoothly with IgNj^, i.e. the discrepancy between 
(b) and (a) rises with increasing primary energy ( Fig.lc ). 

(iv) Within each of (a) and (b) the differences between the (InA) values using different 
combinations of data ( N^, N* etc. ) are systematic and real ( see (i) ). 

(v) In both (a) and (b), the highest estimated (InA) values are for the cases where N* 
is absent - i.e. not used in the analysis. The lowest values of (InA) are where N* has the 
greatest weight in the analysis, i.e. the number of parameters ( 2 in la and 3 in lb ) is a 
minimum. It is an indication that the higher energy muons ( > 2AGeV ) in the observed 
showers are in the greater dehcit compared with model predictions [pl[| . 

(vi) With respect to (v), taking only the two sets ( one in la and one in lb ) which did 
not use N*, it is interesting to note that the former points are close to our lower variant 
( lower dashed curve ) in la and close to our ’best estimate’ ( upper dashed curve ) in lb. 
The differences are, therefore, signihcant. 

The difference in the (InA) values derived using different observables ( Fig.l ) points 
to inadequacies in the models used for the analysis of the experimental data and provides 
an impetus to correct them. However, as remarked already, it is sure that the corrections 
should not be radical, because the difference between the (InA) is not large, typically 
6{InA) ^ 0.2 — 0.4 or, in terms of the ratio of the mean (InA) values from Figures lb 
and la, an increase from 1.12 to 1.37. Below, we show that the mentioned systematic 
difference between the {InA) values can be an indication that the energy distribution 
between the different shower components is slightly different from that in the models: 
specihcally, the actual mean number of electrons W in EAS appears to be higher, that of 
muons slightly lower and that of hadrons W lower still than in the models. 

It is well known from the present models that for the same primary energy the number 


of muons in nuclei-induced showers is higher than in proton-induced ones. On the other 
hand the number of electrons and hadrons in nuclei-induced showers observed in the lower 
half of the atmosphere is lower than in proton ones |^. If one has an opportunity to 
measure the primary energy, by muons, Cherenkov light or another technique, and hnds 
that the shower has a low or a high N^. ( actually the ratio ^ is important ) the 
conclusion will be that this shower is initiated by a proton or light nucleus. On the 
contrary, if one hnds that the shower has a low the conclusion will be the opposite, 
i.e. that the shower is initiated by a heavy nucleus. This is exactly what is observed in 
the showers at sea level. In order to get a better consistency between the (InA) values 
obtained using N^, Nj[ and ( Fig.l ) the model should have higher Ne/Nj[ and smaller 
Nh/Nj[ ratios. This situation could be achieved if the corrections give a reduction of N^, 
a smaller reduction of NjJ' and an increase in N^. 

The same conclusion on the improvement of the model can be drawn from an anal¬ 
ysis of KASCADE event rates §]. Both the muon and hadron trigger rates, observed 
by KASCADE, are lower than expected from the model calculations. This discrepancy 
indicates again that the actual numbers of muons and hadrons in EAS are lower than in 
the models, although the energy region responsible for the trigger rates is lower than that 
analysed for the mass composition around the knee. The argument that the lower trigger 
rates might be due to a lower primary intensity is disproved by the fact that the ratio of 
hadron to muon trigger rates, which is independent of the absolute intensity, is also lower 
than in the calculations [^|. To get agreement the predicted muon trigger rate should be 
reduced by a factor 0.89 ± 0.06, and the hadron rate by 0.54 ± 0.08 compared with the 
rates given by the QGSJET model. This analysis is of value because it indicates that the 
needed reduction of the number of muons in the model should be about 6%, for hadrons it 
is bigger - about 29%. Because muons and hadrons are the products of hadronic cascades, 
it is evident from the energy balance that to reduce the energy contained in the hadronic 
cascade one has to increase the energy transferred to the electromagnetic cascade. At this 
stage it can be remarked that an increase in this energy is predicted by theory ( although 
it appears not to have been included into the models hitherto ). 
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4. Numerical estimates 

4.1. The balance between the EAS components 


To evaluate the effect of the proposed change of the balance between different cascade 
components we applied the semi-quantitative analytical approach. We assumed that 
there are only nucleons (N), pions (vr), muons-|-neutrinos (/ii/) and the electromagnetic 
component (ey) in the cascade. The total and partial inelasticity for nucleons: and 

and for pions = 1 and interaction mean free paths for nucleons A at and 
for pions were taken to be energy independent, i.e. constant. The system of kinetic 
equations, describing the longitudinal development of the nucleon, pion, muon-|-neutrino 
and the electromagnetic energy of the cascade has been taken as 
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Here X is the atmospheric depth, A^r = Att = -^, Ecr - the critical energy for pion 
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decay in the air, {E^^) - the mean energy for pions. The boundary condition at X = 0 is 


Ett + E^^ + Egj — 0 
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where Eq is the total energy of the cascade. The solution of the system of equations is 
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We calculated also the longitudinal development of the electron size of the shower W 


as 


'Ve(-’f) = lo - y)dy 

The longitudinal profile of the electromagnetic cascade has been taken as 
iVe(t) = -^j^exp{t — tmax — 2tlnS), where the atmospheric depth y is measured in units of 
the radiation length to: i = ^, tmax = 1-7 -|- 0.76/n(-^), S = — |Q. Here Eg^ is the 

total energy of the electromagnetic cascade, to = 37.1 gcm“^ and [3 is the critical energy 
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for the electromagnetic cascade in air, equal to 0.081 GeV. Although this approximation 
is recommended for nuclear cascades, its use for the electromagnetic cascades is more 
appropriate for our energy-balance approach, since it gives the correct estimate of the 
depth of the shower maximum ( X^ax ~ 610gcm~^ ) and its elongation rate ( ER 
65gcm~^ ). 

The integrals for E^^, and E^.^ were calculated numerically. The basic parameters 
for the calculation were taken as A'j^j=0.60, A';|^=0.20, ^^^=0.33, A7v=90 gcm“^, A^=120 
gcm“^, Ecr = 90 GeV, {E^,) = 300 GeV. 

The results of the calculation are shown in Figure 2 by the full line. According to 
our suggestion, we increased the energy fraction transferred by nucleons into the electro¬ 
magnetic component, leaving all the other characteristics untouched. We increased 
from 0.20 to 0.26 ( the reason for this particular value will be clear later ). The result 
is shown in Figure 2 by the dashed line: the muon energy at sea level ( 1000 gcm“^ ) 
decreased by ~6%, the hadron energy decreased by ~23%, the energy transferred into 
the electromagnetic component increased by ~2%. 

4.2. The triangle diagrams 

The balance of the energy contained in the major components of the shower is con¬ 
venient to analyse using the so-called “triangle diagrams” If the height of an 

equilateral triangle is equal to 1, then for each point inside this triangle the sum of the 
distances to its sides is equal to 1. If we know the energy fractions carried by the electro¬ 
magnetic (^67)) muon (Sfj_) and hadron (Sh) components of the shower at the observation 
level, so that + 6^ = 1, then each shower can be presented by a single point 

inside the triangle. Our basic shower ( Figure 2, full line ) is shown in Figure 3 by a full 
circle ( 5e7=0.394, <5^=0.505, 5/i=0.101 ). The desired direction for the shift of the energy 
balance in the modihed model is shown by the straight arrow in Figure 3b. 

However, despite the increase of the energy transferred into the electromagnetic com¬ 
ponent, the preserved electromagnetic energy and the electron size of the shower Ng at 
sea level decreased by 14% due to its faster development and then faster attenuation of 
the cascade. The point in the triangle diagram moved in the different direction ( open 
circle ). The shift is similar to the case when we increase the total inelasticity and 
with Kf; = 
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LONGITUDINAL DEVELOPMENT OE 1 PEV CASCADES IN THE ATMOSPHERE 

Development of the nucleon (N), pion (tt), muon (fu.) and electromagnetic components {ey) 



0 200 400 600 _ 800 1000 

Atmospheric Depth, g + cm 



0 200 400 600 _ 800 1000 

Atmospheric Depth, g*cm 


Figure 2: The longitudinal development of 1 PeV cascades in the atmosphere: (a) Fractions of the 
total energy carried by nucleons (N), pions (tt), muons (/i) and transferred into the electromagnetic (ey) 
component ; (b) Electron Size of the shower Ng. Basic parameters of the calculation are given in the 
text. Full line: = 0.20, ER = 65gcm~'^; dashed line: = 0.26, ER = 65gcm~'^; dotted line: 

= 0.26, ER = 71gcm-^. 
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TRIANGLE DIAGRAM FOR THE ENERGY BALANCE OF 1 PEV FAS AT SEA LEVEL 




Figure 3: Triangle diagram for a 1 PeV shower at the sea level: Sh,Sfj., Se-y are fractions of energy carried 
by the hadron, muon and electromagnetic components of the shower, respectively. • - K;^ = 0.20, = 

0.60, ER = 65gcm-2; Q - = 0.26, = 0.66, ER = 65gcm-2; A - = 0.26, = 0.66, ER = 

7\gcm~‘^. (a) Large scale diagram with the inset magnifying the indicated region, (b) Small scale part 
of the full diagram. The straight arrow to the left of the full circle indicates the desired modification of 
the EAS energy balance. 
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If we were to assume a larger fluctuation a{K^) compared with those in the model this 
would result in a wider distributions of N^, and Nh and in a bigger value of a{lnA). 
Within our approach we cannot give the numerical estimates, but we expect that the 
mean energy balance would not move in the desired direction and the mean value {InA) 
would not change significantly. Therefore, an increase of the mean value or increased 
fluctuations of cannot give the required result. 

4.3. Data from mountain laboratories 


We should remark that the conclusion about the decrease of W relates only to mea¬ 


surements at sea level. At mountain altitudes the increase of Kfl results 


m an increase 


in Ng as seen in Figure 2, however it is not enough here, either. The best complex EAS 
array at mountain altitude was at Tien-Shan and the data there also indicated the need 


for the model to be improved |^. The problem was with the so-called ’observed’ mass 
composition, i.e. the mass composition of primary particles responsible for showers at 
the observation level selected by one of their parameters. The Tien-Shan array is at a 


depth of 690 gem 


-2 


I.e. 


PeV-showers are definitely observed beyond their maximum. 


The A'^e-selected showers have to be enriched by showers initiated by protons and light 
nuclei, while in contrast the A^^^-selected showers have to be enriched by heavy nuclei ini¬ 
tiated showers. However, in fact, the ’observed’ composition at Tien-Shan was about the 
same for these two selections [|^, The consistency can be restored by the same way 
as in KASCADE, i.e. by changing the balance between the muon and electromagnetic 
components in the models, in favor of the latter. 


4.4. Slowing down the cascade 

We have concluded that the mere increase of the energy fraction transferred into 
the electromagnetic component does not change the balance properly, although it does 
help in some respects and, as we have stated, it must be present. We now argue that 
another effect that should be present will bring about the desired effect: the slowing down 
of the development of the cascade in its initial stages. The physics behind this modiheation 
will be discussed later. Here we just demonstrate the validity of this argument in a semi- 
qualitative way. For illustration purposes we slow down the development of the hadronic 
and electromagnetic cascade by increasing the elongation rate ER from 65 gcm“^ to 71 
gcm“^, preserving = 0.26. The result is shown in Figure 2 by dotted lines and in 
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Figure 3 by an open triangle. The direction of the shift is correct and its magnitude can 
be adjusted to achieve consistency in the mass composition. 

Thus we conclude that the increase of the energy transferred into the electromagnetic 
component ( ) combined with the slowing down of the development of cascades in 

their initial stages is the most realistic way to improve the particle interaction model and 
to achieve a consistent estimate of the primary mass composition. 

5. Theoretical arguments 

5.1. General remarks 


All the arguments hitherto were purely phenomenological. However, there are theoret¬ 
ical arguments too which lend support to the phenomenological consideration as will be 
demonstrated. At the moment nearly all the models for hadron-nucleus {hA) interactions 
implemented in CORSIKA are well developed and theoretically justified. Nucleus-nucleus 
{AA) interactions are considered within the framework of the Glauber approach, which 
reduces the AA-collision to the sum of a few hA-collisions. No new processes are assumed 
to occur in AA-collisions. However, there are processes which must be present but which 
have not been included hitherto, and we think that both effects: the increase of and 
the slowing down of the cascade development refer to processes, intrinsic for AA-collisions, 
which will cause at least near-equality of the (InA) values in Figures la and lb. 


5.2. The increase of 

At PeV energies, when the Lorenz-factor of the projectile nucleus approaches ~ 10® 
and its charge Z is ~ 10 — 30, the density of virtual photons in the contracted coulomb 
held is so high that when the projectile collides with a nitrogen ( Z = 7 ) nucleus of the air 
and the impact parameter is small ( central collision ) the probability of electron-positron 


pair production including multiple pair production becomes quite appreciable ^ 


Using the numerical estimates made by Erlykin for p-air collisions, we can extrapolate 
them to Fe-air collisions and expect an increase of for the average collision at PeV 

energies by 0.05-0.13. For specihc central collisions the effect will be even higher. Our 


AA 


0.06 of §4.1 was thus rather conservative. We also emphasize the fact that ordinary K, 
value in AA-collisions is not very high. Typically for Fe-air collisions the average number 
of wounded nucleons n,„ is about 7-8 and 


p^PAn^ ~ 0.03. Therefore even a 

1 ^Fe 


few percent increase of in the hrst interaction can be essential for the subsequent 
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development of the atmospheric cascade. 

Experimental results conhrm the intense production of electron-positron pairs in AA 
collisions. In the CERES/NA45 experiment the number of e^e“ pairs in the low invari¬ 
ant mass region rUe+e- < 200MeE/c^ exceeds the number in the higher mass region of 
rUe+e- > 200MeE/c^ by a factor of ~ 10^. For the CERES/NA45 collaboration it is just 
a background, because they are interested in dileptons with high invariant mass ER ^ 


as a signature of an elusive quark-gluon plasma. We believe that to check our assumption 
the study of the dilepton production should be concentrated preferentially in the region of 
invariant mass below lOMeV/c^. 

An additional energy transfer into an electromagnetic component can arise also from 
an excess of direct photons, which has been predicted theoretically as a signature of 


the quark-gluon plasma |^| and is now observed in AA-collisions both at low and high 
transverse momenta |^, 

Another theoretical possibility is to increase the transfer of the energy into the elec¬ 
tromagnetic component, by way of hA-interactions; this could be connected with the 
breaking of the diquark and its subsequent recombination into 3 7r°-mesons P6| , ^ . 

5.3. Slowing down of the cascade development 

The theoretical motivation of the second assumption is also connected with the prop¬ 
erties of AA-interactions, but relates mainly to peripheral collisions. It is known that even 
at high energies a projectile nucleus does not fully disintegrate into constituent nucleons 
in the hrst interaction but fragment into few pieces of different mass [^. Some of these 
fragments are excited and, after the de-excitation, if it occurs in space, give rise to MeV 


gamma-ray lines, observed from ’discrete’ sources and the interstellar medium [^. The 
lifetime of the excited fragments varies from a ’nuclear’ time ~ 10“^^ sec to millions of 
years. For AA-interaction at PeV energies both the lifetime before de-excitation and the 
energy of emitted gamma-quanta are extended by the factor of 10® — 10® due to rela¬ 
tivistic effects. As a consequence one can expect an additional sub-PeV electromagnetic 
cascade to be initiated a few hundred meters below the point of the hrst interaction. 
This effect will slow down the development of the electromagnetic cascade and shift its 
maximum. 

In general the shift is limited to no more than 1-2 mean free paths for the inelastic 
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interaction of the fragments. However, there might be an additional effect which is, it 


must be admitted, speculative. It is based on the so-called sling” effect |^, |^. and is 
connected with the deformation of the shape of the nuclear fragment, which rotates with 


a high spin . As a result of such deformation one can expect an increase in the mean 
free path for inelastic interaction which enables the excited nuclear fragments to penetrate 
deeper into the atmosphere, slowing down the cascade development and resulting in a 
further shift of X^ax- There will be an additional increase of N^. for nucleus-induced 
cascades like that shown in Figure 2b by the dotted line. 

All the consequences of the ’sling’ effect: polarisation of the secondary nuclear frag¬ 
ments, deformed shape of the nuclei rotating with a high spin and cross-section fluctua¬ 


tions of the excited nuclei are true effects observed at low energies |^. The problem is 
whether they still hold at higher energies. However, even if they don’t, the de-exitation 
of the nuclear fragments itself will give the needed shift of X^ax if their lifetime is not 
less than 10“^^ sec. 

It is difficult to make reliable numerical estimates of the shift since there is not enough 
experimental data and no good theory of the effect. We think that to check our assumption 
the study of the properties of the secondary nuclear fragments should be made in high 
energy AA-collisions with different nuclear beams. 

6. Discussion 

6.1. Application of the new model to PeV energies 


6.1.1. The mass composition problem 

In order to get better consistency of the results on the primary cosmic ray mass 
composition, the combination of two effects is needed. The increase of only is not 
enough because although it results in a decrease of and Nh it also leads to a fall in Ne if 
the measurements are made at sea level, whereas we need an increase in Ng. The slowing 
down of the electromagnetic cascade alone does not change the ratio between muons and 
hadrons. The slowing down of just the nuclear and electromagnetic cascades ( e.g. by 
increasing the mean free path for nucleon interaction Aat without increasing ) results 
in an increase in both Ne and Nh. Only the combination of both effects: the increase of 
and the slowing down of the development of the electromagnetic cascade changes the 
energy balance in the needed direction: and Nh decreased, Ng and X^ax increased 
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( Fig.3 ). The required changes of and Xmax are small, because the sensitivity of 
the energy balance to these changes is rather high. 

The result should be a coming together of the {InA) values in Figures la and lb. 
There will be an increase in the values in Figure la and a small reduction in those in 
Figure lb. In view of the problems associated with N* ( see §3 ) most of the predictions 
in Figure lb should be disregarded; the ’best estimate’ is then just a little below the 
crosses in Figure lb - these values are reassuringly close to our latest estimate indicated 
by the upper dashed curve in Figure 1. 

6.1.2. The Xmax controversy 

Although the Xmax aspect for EAS is one of some complexity we discuss it briefly. 
The assumption of a slowing down of the cascade development and the consequent shift 
of Xmax in AA-interactions helps us to understand the striking difference between the 
mass composition derived recently by Swordy and Kieda p and Fowler et ah [§], using 
Cherenkov light measurements and those by Roth et ah 0 and others, without. Those 
with Cherenkov light give a significantly lighter mass composition and lower intensity of 
the primary energy spectrum than those obtained without the use of Xmax 0]- If there 
is indeed a shift of Xmax and a consequent increase of Nf, the shower looks more like that 
initiated by proton and its primary energy, based on the on-ground measurements is un¬ 
derestimated. Insofar as the Nf. values are most in doubt, from the prediction standpoint 
the real mass composition should be closer to that derived by using direct measurements 
in the stratosphere or using EAS muons and hadrons on the ground ( Fig. lb ). 

6.2. Application of the new model to higher energies 

If the increase of and the slowing down of the development of the cascade are 
due to electron-positron pair production, and de-excitation of nuclear fragments, there 
will be signihcant further consequencies of the model because of its charge and energy 
dependent QED effect and the displacement of the depth of shower maximum. Both these 
effects will grow with increasing collision energy. In this case it is not surprising that the 
inconsistency in the estimates of the primary mass composition grows with energy, as 
noticed by Roth et ah [0 and indicated in Figure Ic. 

Of greatest importance is the mass composition at the highest energies, where, due to 
the low anisotropy and the rapid growth of Xmax with energy, particles are generally con- 
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sidered to be extragalactic (EG). There is controversy, here, between those who subscribe 
to the common view that the EG particles are protons ( eg [^, |^, ^ and ourselves ( eg 
following much earlier work by Tkaczyk et al. [^). A later publication will deal 


with this aspect in detail but a few brief remarks here are in order. 

(i) The change of will be even bigger at the highest energies ( say above 
GeV ) and the result will be an increase in the mean primary mass: this follows from the 
analysis of Ave et.al.[^|, in which inclined muon showers were detected - the correspond¬ 
ing primary energy will now be higher and the rates more in accord with expectation 
for a signihcant fraction of primary iron nuclei, a recent recalibration of Haverah Park 
energies 0 leads to a reduction in primary energy and an even further increase in the 
iron fraction. 

(ii) The displacement of the depth of maximum to greater values will, again, lead to an 
increase in the mean primary mass as indicated by the X^ax values from the experimental 
data ( see Fig. lb ). 

7. Conclusions 


Our analysis of the present situation with the high energy interaction models indicates 
that there is no support for the introduction of a radical change of the models, sometimes 
proposed to explain the origin of the knee in the primary cosmic ray energy spectrum. 
However, the inconsistencies in the interpretation of the experimental data on the primary 
mass composition, obtained when different EAS components are used for the analysis, 
indicate the need for some improvement to the models which were used hitherto. We 
propose that the most promising way is to introduce an additional ( a few percent ) energy 
transfer into the EAS electromagnetic component combined with a slowing down of the 
cascade development at its initial stages, which is followed by a small ( ~ 20 — 30gcm~‘^ ) 
shift of Xmax into the deeper atmosphere and the consequent increase of Ng. The most 
likely processes which can be responsible for such changes ( e.g. electron-positron pair 
production including multiple pairs, direct photons from the hypothetical quark-gluon 
plasma, excitation of the secondary nuclear fragments ) are those which occur in nucleus- 
nucleus collisions and they should indeed be present at some level. The importance of 
these processes is expected to grow with energy and offer the hope of resolving some 
controversies at very high energies. 
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